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Abstract 

The theory of photon emission, absorption and scattering In a 
relativistic plasma of positrons, electrons and photons Is studied. 
Expressions for the emlsslvltles and absorption coefficients of pair 
annihilation, pair production and Compton scattering are given and evaluated 
numerically. The conditions for negative absorption are Investigated. In a 
system of photons and e + -e" pairs, an emission line at ~ 0.43 MeV can be 
produced by grasar action provided that the pair chemical potential exceeds ~ 
1 MeV. At a temperature of ~ 10% this requires a pair density > lO^cm"^, a 
value much larger than the thermodynamic equilibrium pair density at this 
temperature. This emission line could account without a gravitational 
redshift for the observed lines at this energy from gamma ray bursts. 
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I. Introduction 

The 'mission and absorption of photons In cosmic sources are governed by 
many processes. At temperatures of the order 10® to 10^K, typical of gamma 
ray burst sources, two of the most Important ones are pair production and 
annihilation (e + +e" ♦ y + y) and Compton and Inverse Compton 
scattering (e+Y ♦ e'+Y'). Arguments based on the observed photon Intensities 
of gamma ray bursts and the likely distances and sizes of their sources, lead 
to the conclusion that the source regions of at least some of the bursts are 
optically thick (Cavallo and Rees 1978, Schmidt 1978). 

Photon absorption In yy pair production has been discussed In the 
literature (Gould and Schreder 1967), but no calculation has Included the 
effects of the stimulation of the annihilation or the suppression of pair 
production due to large photon or particle occupation numbers. When these 
stimulation and suppression effects are taken Into account, the possibility 
exists for negative absorption (Varma 1977). The condition for this Is a 
population Inversion, which In the present context Is a pair density that 
exceeds the thermodynamic equilibrium density. 

A recent review of gamma ray burst observations has been given by Cline 
(1981). Of particular interest for the present paper Is the existence of an 
emission line seen from several ganmta ray bursts In the energy range from 0.40 
to 0.46 MeV (Mazets et al. 1979, Teegarden and Cline 1980, Mazets et al. 

1981). These lines are probably due to e + -e” annihilation radiation. If so, 
e + -e - pairs should be present In large numbers In the burst sources, and the 
sources should be sufficiently hot to produce the pairs, but the source 
regions should not be In thermodynamical equilibrium because no lines can then 
be seen. We aim the calculations of the present paper to astrophyslcal sites 
where such conditions might exist. 
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We provide the basic analytic treatment In Section II, we give numerical 
results In Section III, we discuss some astrophyslcal applications In Section 
IV, and give our conclusions In Section V. 

II. Emlsslvltles and Absorption Coefficients In a Relativistic Plasma 
We consider systems characterized by temperatures of the order of the 
electron rest mass energy In which photon-photon collisions can produce much 
larger pair densities than the ambient electron densities of the astrophyslcal 
sites of Interest. We therefore consider only cases In which the electrons 
and positrons have equal densities. As convenient analytical expressions, 
which allow both equilibrium and non-equilibrium situations, we use Bose- 
Elnsteln distributions for the photons and Fermi -Dirac distributions for the 
pairs (see Landau and Llfshltz 1958). The reference frame In which these 
distributions are Isotropic Is designated as the plasma frame. 

We assume equal temperatures for the positrons and electrons, T + = T_ = 

T t , but allow the photon temperature, T y , to differ from T t . Since the e + and 

e" densities are equal, n + * n_ * n ± , these particles must also have equal 
chemical potentials, u + = u_ = u t . The photon chemical potential, u Y , is zero 
for a blackbody distribution. We allow non-blackbody photon distributions, 
but only zero or negative values may be assigned to Uy The pair chemical 

potential can be positive, zero or negative. If u ± * 0 the pairs are In 

thermodynamic equilibrium with blackbody photons. 

In terms of these temperatures and chemical potentials, the photon and 
pair densities can be written as (e.g. Landau and Llfshltz 1958) 


n - (AwVf 1 /^ n , 


( 1 ) 
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where p Is the photon or particle momentum and n Is the occupation number. 
These are given by 

n y • {expL (E y -m y )/ kT y J-l ) -1 (2) 

for the photons, and by 

n ± = {expl(E ± -u ± )/kT ± J + l} -1 (3) 

for the particles, where E y Is the photon energy and E t Is the particle total 
energy (kinetic plus rest mass). Equation (1) counts both polarization states 
for the photons and both spin states for the particles. 

The blackbody (y Y = 0) photon density Is shown by the line y In 
Figure 1. This quantity has the simple analytic form iiy = 2c(3)w" 2 (kT/c1 l i) 3 , 
where c(3) * 1.2021.... The other curves In this figure show pair densities 
for various chemical potentials u±. Pair densities corresponding to positive 
or negative u t yield higher or lower pair annihilation rates, respectively, 
than the pair production rate of blackbody radiation. The Fermi -Dirac 
distributions (equation 3) tend to Maxwel 1 -Bol tzmann distributions in the 
limit of large -u t . 

We proceed to define the photon emlsslvlty and absorption coefficient for 
pair production and annihilation. In particular, we are Interested In 
obtaining a correct expression for stimulated annihilation which has not been 
taken Into account In previous treatments of absorption In photon-photon pair 
production. Stimulated emission has, of course, been taken Into account for 
other processes (e.g. Bekefl 1966). But we cannot use the standard 
expressions for emlssltlvles and absorption coefficients because pair 
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production and annihilation dc not fit the usual pattern In which photons are 
emitted or absorbed singly and the matter has the same form before and after* 
events. We therefore proceed as follows: 

For two-photon annihilation and pair production, e + +e" ♦ -rj+Y 2 » the 

transition rate In vacuum In either direction between the photon states In 

3 3 3 3 

d Pjd pg and the pair states In d p + d p_ can be written as 

w = (4wV)" i 2 d 3 p 1 d 3 p 2 d 3 p + d 3 p_ 6 4 (p 1 +p 2 -p + -pjx. (4) 

Here j>j and p 2 are photon momenta, p + and p_ are momenta of the pair, pj, p 2 , 
p + and p_ are the corresponding 4-momenta, and X Is proportional to the 
squared matrix element of the Interaction, summed and averaged over spins and 
polarization. 

To obtain the total annihilation and pair production rates we must 
multiply equation (4) by an appropriate combination of the occupation 
numbers. This can be obtained from energy conservation ( Ej+E 2 * E + +E_) and 
the equilibrium condition (uj a and T t * T y ) , where Ej and E 2 are the 
energies of the two photons. Equations (2) and (3) then yield 


n l n 2 ( 1 ” n + ) ( 1 -n_ ) * n + n_(l + nj)(l + n 2 j 


(5) 


Multiplying both sides of this equation by the transition rave (equation 4) 
and Integrating over all four momenta, we obtain 


1 n.t^p n d ^P A- 

7 /.-j~ :" x' i / ^1 

4« Ti 4* 


ann 


i n,o p. Hod Pp m - 

’ T I I -^T 1 “.'l-n.lll-n.lllc 35]^. 


(b) 
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The left-hand side of equation (6) Is the total pair annihilation rate, 
while the right-hand side Is the total pair production rate. The Invariant 
product of the flux factor I and differential cross section do/ do (Jauch and 
Rorllch 1955) Is obtained from X by Integration over all final state variables 
except the angles of one particle: 

[1° cRi^ann = /Pi ^Pi/P2 ^Pl + P2”P+”P-^ 

LIc 7BT ] pp * / P+ 2 d P+ /P_ 2 d P_J d£J _ < Pi + P2 ”P+"P- )x • (8) 

The factor ^ In equation (6) Is Introduced so that each distinct pair of 
photons In either Initial or final state Is Included just once. Note that the 
annihilation Is stimulated by the presence of the bath of photons, while pair 
production Is suppressed by the presence of the bath of pairs. However, 

If = p^<<0, the pair and photon distributions are reduced to 
Maxwell -Boltzmann distributions and this removes all the stimulation and 
degeneracy effects. 

To define a photon emlsslvlty and absorption coefficient It Is necessary 
to Investigate the balance of reactions Involving photons In an Increment 
d 3 pj. We proceed as In the derivation of equation (6), except that the 
Integration extends only over d^, d 3 p + , and d 3 p_. In this case, however, It 
Is no longer possible to express the left side In terms of the annihilation 
cross section, because the necessary Integration of equation (7) has not been 
performed. But we can Interchange the order of Integration and use equation 
(8) Instead. The required balance Is then given by 
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(l+n 1 )d 3 p 1 d 3 p 2 

B ^d 3 Pl ^ d 3 p 2 
4w 3 fl 3 1 4* 3 H 3 


/ <to + n + n.(l+n 2 ) Lie ^yj pp ■ 
/da + n 2 (l-n + )(l-nJLIc^J pp> 


(9) 


Collecting terms In n^, we obtain 


dE 1 dB lV E l ) 


cn^Ej dEjdfl^ 
3 3 — 

4* wer 


W’ 


(1U) 


where the rate of spontaneous emission, J (Ej), and coefficient of linear 
absorption, K^lEj), are given by 


2 3 


PP 


(ID 


and 


K (E.) 
YY 1 


^3 

/ — 7^/ dfi + Ln 2 (l-n + )(l-n_) -n + n_(l+n 2 )][I-g^J . 

4*'Tr pp 


( 12 ) 


In equation (11) the annihilation emlsslvlty Is expressed In terms of the pair 
production cross section unlike the approach where this Is done In terms of 
the annihilation cross section (e.g. Ramaty and Meszdros 1981). In the 
expression for the absorption coefficient (equation 12), the first term In the 
brackets Is due to absorption by the photon bath while the second term Is the 
contribution of Induced annihilation. When the equilibrium condition (5) Is 
satisfied and u Y « 0, equations (11) and (12) yield Klrchhoff's law. 


J /K 

tY YY 


4i 3 (Kc) 3 exptEj/kTM 


i I 


BB 


»E,), 


(13) 



where I BB Is the blackbody Intensity. 

Equations (11) and (12) are valid also for nonequl librium situations 
provided proper nonequl llbrl inn occupation numbers are used. In terms of 
equations (2) and (3), the most general nonequilibrium distributions are 
obtained If T and y^*y ± With such distributions, the total 
annihilation and pair production rates are not equal and, moreover, can 
become negative. While for * T Y and y t * y.y, Is alwqys positive, for 
nonequilibrium conditions the contribution of stimulated annihilation can 
exceed that of absorption provided an appropriate population Inversion takes 
place. By substituting equations (2) and (3) Into equation (12) we find that 
for T ± * T y such an Inversion occurs If 2y^ > y^. In this case i: 
negative for Ej < 2y^-y^ . 

For the system to exhibit grasar action, however. It Is necessary that 
the total absorption coefficient be negative. For the system of photons and 
pairs that we consider here, the only Important process other than pair 
production and annihilation Is Compton scattering. We Ignore the weaker 
processes of bremsstrahlung and double Compton scattering. We note, however, 
that synchrotron radiation could potentially be very Important, but because we 
are free to choose an arbitrarily low magnetic field Intensity, we Ignore 
synchrotron absorption In the present discussion. 

For Compton and Inverse Compton scatterings (y^e ♦ y 2 + *') we proceed In 
essentially the same w«y as for pair production and annihilation. Using p^ 
and p 2 for photon momenta and p and p' for electron and positron momenta, we 
find the equation for the overall balance between Inverse Compton scattering 
and direct Compton scattering to be 
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Ha p Hi J- 


(14) 


where the factor of 2 takes Into account the contributions of both electrons 
and positrons. Here again the Invariant product of flux factor *nd 
differential cross section Is given by Integrating the appropriate squared 
matrix element over all final state variables except the photon angles. For 
the direct scattering we have 

Lie -g£] c - /p 2 2 dp 2 /p ,2 dp7da , « 4 (p 1 +p-p 2 -p*)X c . (15) 

For the Inverse scattering, designated by subscript C', the same form applies 
with the substitutions p 2 ♦ Pj and p' ♦ p. Equation (14) Is valid for any 
choice of jij and \iy as long as T A * T y . 

To obtain expressions for the emlsslvlty and absorption coefficient, we 
again choose an Increment d^pj and Integrate over the other three momenta. As 
for pair production and annihilation, we must Interchange the order of 
Integrations on one side In order to Introduce a cross section. After the 
appropriate manipulations we obtain the Compton emlsslvlty and absorption 
coefficient In the presence of the bath of photons and pairs 


w 


4«^(Hc) 3 




(16) 
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K C <E 1 ) - Z ^“^jW*2 Ln * tl+n z) ( l" n t ) * n 2 , *i a " n t ,]tI i J C * (17) 

The Compton emlsslvlty, equation (16), represents the scatterings of photons 2 
Into the element dE^do^. For the absorption coefficient, equation (17), the 
first term In the brackets Is due to scattering of photons out of dEjdQi, 
while the second term represents the stimulated scatterings of photons 2 Into 
dEjdOj. Kc(E}) can become negative and a necessary condition for this Is T Y > 
Tj. In our subsequent analysis, however, we shall or*y consider system with 
for which K c Is always positive. In such cases equations (16) and 
(17) satisfy a modified Klrchhoff's law 

cE i 2 -i 

J-IE.l/IUE.) ■ , 1 , {«pl(E.-u v )AT 1-1) *. (18) 

For the numerical evaluations shown below we Stave used the expressions of 
Jauch and Rohrllch (1955) for the flux factors and differential cross 
sections. We must also express all quantities In term of Independent 
variables of Integration. For pair production from photons of energies Ej and 
E 2 , the flux factor Is given by 


! pp " U C0S8 12‘ 


(19) 


where 8^ Is the angle between th? two photons In the plasma frame. In the 
center-of-momentum frame (CM frame), the positron (or electron) has Lorentz 
factor y given by 
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Wmc 2 ) 2 -VfeEjEjd-cose^) . (20) 

The Lorentz factor y^m for the tranformatlon between plasma frame and CM 
frame Is 


t CM ‘ <Ei*E 2 ) /( 2yi»c 2 > 


( 21 ) 


T he velocity associ»tcu with ycm Is In the CM fr.sme the direction of the 
positron Is given by colatitude 9 fro.n the direction of photon 1 and azimuth + 
from the plane of b cm and photon 1. The electron direction Is diametrically 
opposite. The differential cross section Is 


do- 

m- 


r c e l-B' + cos 4 e+2Y‘ Z B 2 s1n 2 0 

77 (l-oWe? • 


( 22 ) 


where r Q Is the classical electron radius (£.818x10" l,} cm)and B * (1-y” 2 )^. 
The angle 8j between photon 1 and Bq^ Is given by 


cose^ 



[EjlYQjjjYmc 2 )’ 1 -!], 


(231 


and the angle 6+ between the positron and Bqh Is 


cose + * cosfljcose+slne^slne cos*. (24) 

Finally the energies of the pair In the plasma frame are given py 

E i “ Y CM T * 2(lt8 CM aco#e + , » 


(25) 
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and the occupation numbers are given by equations (2) and (3). Thus 
everything needed to evaluate equations (11) and (12) has been found from E^, 
E 2 , e 12 , 9, and t . A Monte-Carlo technique was used to evaluate the four- 
dimensional integrals over E 2 , e^ 2 , e, and $ . 

For Compton scattering all quantities can be evaluated directly in the 
plasma frame, thanks to expressions given by Jauch and Rohrlich (1955). For 
an initial pnoton of energy Ej and an initial electron or positron of energy 

E = y ± mc 2 = mc 2 (l-B ± 2 f (26) 

the flux factor is given by 

Iq = l-e + roso 1 , (27) 

where is the angle between electron and photon. The direction of the 
scattered photon (photon 2) is given by colatitude 0 from the direction of 
photon 1 and azimuth $ from the plane of photon 1 and the initial electron. 
The angle o 2 between photon 2 and the initial electron is given by 

cos = cos <*j cose + sine^sine cos$. (28) 

The energy E 2 of photon 2 is given by 
E 9 1-e. coso, 

2 ± t (29) 

l+(E 1 /E + Hl-COS0)-9 + COSo 2 

The differential cross section is 
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kj * EjY+U-BjCOSoj ) (31) 

kg * E 2 t ± (1-B ± co^). (32) 


The energy of the final electron Is 


E' « E ± + E r E 2 , (33) 

and the occupation numbers can be found from equations (2) and (3). Again 
everything needed to evaluate equations (16) and (17) has been found In terms 
of Ej, E t , aj, e, and + . A Monte-Carlo technique was used to evaluate the 
four-dimensional Integrals over E ± , a lt 8, and +. 

III. Numerical Results 

We have evaluated equations (11), (12), (16) and (17) for various choices 
of T t , t ^ and y As already Indicated, we limit our discussion here to 
cases with equal pair and photon temperatures, T* * T y s T. We allow, 
however, arbitrary values for uj and y 

We consider first the case of thermodynamic equilibrium, u* * w Y * 0. 

The emlsslvltles and absorption coefficients for this case and T = 3x10^ are 
shown In Figure 2, where * j YY + j c and K t * K YY + Kq. As can be seen, the 
absorption coefficients are positive at all photon energies and Klrchhoff's 
law Is satisfied by all processes. 

We next consider a case of equilibrium between pair annihilation and pair 
production by non blackbody photons. Numerical results for n* * ^ * -2 MeV 
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and T - 3x10^ are shown In Figure 3. Because for these parameters both the 
particle and photon occupation numbers are very small (<10” 3 ), the results of 
Figure 3 closely approximate emlsslvltles and absorption coefficients 
appropriate for Maxwell -Boltzmann distributions. Indeed, the emlsslvlty jyy 
shown In Figure 3 Is essentially Identical with the annihilation emlsslvlty 
calculated by Ramaty and Mdszdros (1981) using a Maxwell -Boltzmann 
distribution and the pair annihilation cross section. As discussed by these 
authors, as well as Zdzlarskl (1981) and Aharonlan, Atoyan and Sunyaev (1980), 
the peak of the annihilation emlsslvlty occurs at a higher energy than mc^ ■ 
0.511 MeV, because the annihilation photons must carry away the kinetic 
energies of the pairs In addition to their rest mass energy. This effect Is 
very obvious In both Figures 2 and 3. 

From the emlsslvlty jyy of Figure 3 we can evaluate the full width at 
half maximum (FWHM) of an optically thin annihilation feature produced In a 
Maxwellian plasma of Tj, * 3 x 10 9 K. We find that for this temperature the 
FWHM (~ 850 keV) Is In excellent agreement with the calculations of Ramaty and 
Meszaros (1981) who deduced the dependence of FWHM on T t for such plasmas. 

The absorption coefficients of Figure 3 are positive at all energies. While 
Klrchho*'' taw Is not satisfied (because u y *0), the J's and K's of Figure 3 
do satisfy the modified Klrchhoff's law (equation 18) for all processes. 

We turn now to the study of cases with Inverted populations. For T^ * T y 
the Inversion condition Is > u^/2. Since <0, the Inversion threshold 

Is at a value of u* that Is higher than or equal to the value at which pair 
production and annihilation are In equilibrium, * u y . The two values are 
equal for blackbody photons, * 0. When an Inversion occurs, Kyy Is 
negative for Ey < 2^-Uy, but since T^ ■ Ty, K c Is positive for all Ey. 

Grasar action can occur only If ■ Kyy + K^<0. Since Kq Is proportional to 
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n ± while the portion of K n due to stimulated annihilation varies as n^ 2 , a 
sufficiently large density Is needed for -K^ to exceed K c . This Implies a 
threshold for u t which Is higher than the threshold required for just K yY to 
be negative. 

To Investigate this threshold we have evaluated K t * K YY + K c as a 
function of for given temperatures and u^. We have carried out 

Q 

calculations In the temperature range 0 < T < 5x10 K, where the lower limit 
corresponds to fully degenerate electrons and positrons. We flr.d that the 
threshold for grasar action Is close to * 1 MeV and does not depend 
strongly on u Y and T. As can be seen from Figure 1, this corresponds to a 
pair density threshold of a few times 10 30 cm - ^. 

We show numerical results In Figure 4 for T * 3xlO^K, = 0 and * 1.1 
MeV. As can be seen, K yy Is negative for E y < 2u t * 2.2 MeV and positive at 
higher energies. K c Is positive at all energies and Klrchhoff's law (j c /K c * 
I BB ) Is satisfied for Compton scattering since = 0. In the energy range 
from about 0.25 MeV to 0.7 MeV, -K n exceeds K c and hence K t Is negative. If 
the source Is optically thick and K t Is negative over a sufficiently large 
spatial region, then the radiation Intensity has a sharp peak at a photon 
energy at which -K t Is maximal. The value of this energy, - 0.43 MeV from the 
numerical calculations, Is determined primarily from the energy at which -K^ 

Is maximum, shifted somewhat according to the slope of K c at that energy. 

From Equation (12) we can also express In the form 

3 

vv * • w * ^^ da * , '2 ll - n . ,n -''- ltI ^ ] pp • <34) 


Above the grasar action threshold, the first term, due to stimulated emission. 



17 


Is much larger In magnitude than the second term which 1$ due to absorption. 
As discussed above, J YY (E Y ) Is broadly peaked at an energy greater than me 2 , 
reflecting the kinetic energy of the annihilating pairs. The division by E y 2 
(from the factor of density of states) shifts the peak to an energy somewhat 
less than me 2 . 

The most extreme case of population Inversion arises when T Y * T ± = 0. 
The system then Includes no thermal photons at all. The particle states are 
fully occupied or degenerate (n^ = 1) up to the Fermi momentum pp, and are 
empty ( = 0) above pp. The Fermi momentum Is related to the chemical 
potential by 

cp F = (u ± 2 - m 2 c 4 ) 1/2 , (35) 

and to the pair density by 

n ± = OwVj^pp 3 . (36) 


At T y = 0 there Is no absorption term In K yY because there are no photons 
In the bath. Therefore K YY can only be zero or negative and Is proportional 
to j Yy : 


V, ,s -7r^W' 


(37) 


Both J YY and K YY are non-zero only between the kinematic limits 


u ± - cp F < * epp . (39) 

These limits correspond to the annihilation of a pair with equal momenta of 
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magnitude p F Into one photon with momentum parallel to the pair momenta (upper 
kinematic limit) and another photon with momentum antiparallel to the pair 
momenta (lower kinematic limit). 

Similarly, j c * 0 and there Is no stimulated emission term In Kg because 
there are no photons present to be scattered to contribute these effects. 

Only the absorption term In Is present, so that Kq Is necessarily 
positive. For small E y , only those particles with momentum near Pp can 
contribute to K c , because the particle can only be scattered Into a previously 
unoccupied state. 

Numerical results are shown In Figure 5 for T * 0, = 0 and 

* 0.85 MeV. Here the broadening of Is caused by the motion of the 

degenerate particles even though their temperature T t s 0. We shall refer to 

this effect as degeneracy broadening. By evaluating for other values of 

as well, we find that the FWHM for degeneracy broadening Is proportional to 

1/3 

the Fermi momentum pp, and hence proportional to n ' for both non- 
relatlvlstlc (p F « me) and relativistic (pp > me) degenerate distributions. 

It Is given by 

FWHM - (4.3 x 10" 8 keV) n ± 1/3 (40) 

when n t Is expressed In cm -3 . 

In Figure 5 the klnematlcal limits on j Y and K yy are quite evident, as 
Is the suppression of Kq by degeneracy at low photon energies. The peak of 
J n Is at about 0.75 MeV, but because of the E y " 2 factor In equation (34), the 
peak of -Kyy Is shifted to about 0.50 MeV. Because of the steep positive 
slope of Kq, the peak of In Figure 5 Is further shifted to about 0.42 MeV. 

From the numerical calculations for other values of T and 
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to<T< bxl0 9 K, 0.8<w t < 1.2 MeV) we find that the peak of -K t Is Invariably In 
the energy range from about 0.40 to 0.47 MeV. When Is close to the 
threshold for grasar action, the peak energy of -K^ decreases with Increasing 
temperature, but the varying slope of Kq resulting from the degeneracy tends 
to cancel this effect. Thus, there Is not much variation of the peak energy 
of -K t with either T or u ± . For the parameters considered, we have found no 
case when the peak energy of -K t equaled or exceeded 0.5 MeV. 

III. Astrophyslcal Applications 

Gamma ray burst sources are likely astrophyslcal sites where the results 
of the present paper could be applied. The large photon densities expected 
(Cavallo and Rees 1978, Schmidt 1978) In these sources should lead to high 
pair production and Compton opacities. The observation (Mazets et al. 1979, 
Teegarden and Cline 1980, Mazets et al. 1981) of an emission line In the 
energy range 0.4 to 0.46 MeV, believed to be due to e + -e" annihilation 
radiation, Is evidence that e + -e" pairs do Indeed play an Important role In 
the physics of gamma ray bursts. But It Is not Immediately obvious how a 
relatively narrow emission line Is produced In a hot and optically thick 
source region. 

Ramaty et al. (1980) and Ramaty, Ungenfelter and Bussard (1981) have 
studied this problem and discussed an optically thin model for the transient 
of March 5, 1979 (Barat et al. 1979, Cline et al. 1980, Evans et al. 1980) 
from which an emission line was observed (Mazets et al., 1979) at ~ 0.43 
MeV. In this model, the ~ 0.43 MeV line Is formed at an energy > 0.511 MeV In 
the last optical depth of the source region by the annihilation of e + -e“ pairs 
that have been cooled by synchrotron radiation prior to their annihilation. 

The shift from above mc^ to the observed energy 1$ due to the gravitational 
redshlft of a neutron star. The observed upper limit on the 
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width (FWHM < 0.2 MeV, Mazets et al. 1979, 1981) Implies a temperature less 
than 3x10% (Ramaty and Meszaros 1981). The upper limit on the width also 
Implies an upper limit on the density because of degeneracy broadening. Using 
equation (40) we obtain n ± < 7xl0 2 ®cm" 3 for the e + -e" annihilation region If 
the line Is produced by annihilation In the last optical depth. 

The density n ± can also be directly calculated (Ramaty et al. 1981) from 
the observed line fluence ($ » 10 photons cm" 2 , Mazets et al. 1979). Let 
R/(n ± ) 2 » 7.5xl0“ 15 cm 3 sec" 1 be the annihilation rate coefficient at 3x10% 
(Ramaty and Mdszdros 1981), A the area of the emitting region, at the time 
Interval In which the observed fluence Is produced and d * 55kpc the distance 
to the source. Then If the line Is formed In a layer of unit optical depth to 
Compton scattering, 

* « 2R Kj.' 1 Aat (4*d 2 ) -1 (41) 

Since R varies as n ± 2 and K c -1 as n ± -1 , ♦ Is proportional to n^. With the 
above numerical values, Kq from Figure 3, and n ± < 7xl0 2 ®cm" 3 , Aat should 
exceed l.SxlO^n^sec. This condition Is well satisfied If the annihilation 
line Is produced over the entire surface of a neutron star, A * 10* 3 cm 2 , and 
during the entire Impulsive phase of the March 5 event, At * 0.15sec (Cline et 
al. 1980). But the optically thin model would face considerable difficulties 
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If future measurements should Indicate that the 11:<e Is narrower than 0.2 MeV, 
or If Aat should turn out, for other reasons, to be small er than 1.5xl0 9 cm*sec 
(e.g. If the emitting area covers only a fraction of the polar cap area). 

The advantage of producing the annihilation line by grasar action Is that 
a narrow line can form In a hot optically thick region. To Illustrate this, 
we have evaluated the photon Intensity perpendicular to a slab of thickness L 
In which the emlsslvlty and absorption coefficient do not depend on position: 

I - (j t /K t ){l - exp[-K t L]} (42) 

Using the j t 's and K t 's of Figure 4 (T + * T y * 3xl0 9 K, * 0, y+ * 1.1 MeV), 

we show In Figure 6 the dependence of I on photon energy E^ and slab thickness 
L. As can be seen, grasar action can Indeed narrow the line. For example, If 
L>10" 5 cm, the width Is less than 0.2 MeV. In comparison, the thermal width Is 
~ 0.8 MeV In an optically thin Maxwell -Boltzmann gas at 3xlo\ (Ramaty and 
Meszdros 1981, or see In Figure 3). 

The peak energy of the annihilation line formed by grasar action Is In 
the range 0.4 to 0.47 MeV, l.e. close to the observed peak energies. Thus, 
the gravitational redshlft of the line due to the compact object which 
presumably produces the burst should be quite low, z < 0.1. This Implies that 
gamma ray burst sources with observed e*-e~ emission lines could be objects 
other than neutron stars, or If they are neutron stars, these stars should 
have small masses (M < 0.6M o , Borner and Cohen 1973). 

Returning to the results of Figure 6, we note that the photon Intensities 

at the centers of the narrow lines entail very large photon occupation 

3 3 2 -2 

numbers. At 0.43 MeV, these numbers, given by n^ » 4n Ire UE^), are 
■ 200 and 8000 for L * iO“ 5 cm and 2xl0"®cm. The development of such high 
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photon densities clearly should modify the j's and K's used to calculate the 
Intensity of Figure 6, but we defer the Investigation of such a nonlinear 
system to a subsequent study* 

Nevertheless, assuming that an Inverted layer with parameters as In 
Figure 6 did exist In the March 5 burst source region, the line fluence + can 
be calculated as follows 

4 - AAtd " 2 /I(E Y )dF T (43) 

Using the results of Figure 6 with L * 10“^cm and 4 <* 10 photons cm" 2 , 
equation (43) yields Aat * 1.2 x 10 6 cm^sec. By comparing with the minimum 
Aat deduced above for radiation produced In the last optical depth (Aat > 1.5 
x 10 9 cm 2 sec), we see that not only can grasar action produce a narrow line 
In a much hotter region, but that the observed line Intensity and width are 
consistent with a much smaller source and/or a much shorter line formation 
time. 

V Conclusions 

We have carried out a fully relativistic treatment of pair production and 
annihilation and Compton and Inverse Compton scattering In a medium containing 
photons, positrons, and electrons, with equal e + and e~ densities. In the 
calculation of the emlsslvltles and absorption coefficients we have Included 
the stimulation of transitions caused by the Bose-Elnsteln nature of the 
photons and the suppression of transitions due to electron and positron 
degeneracy. We have shown that for systems In thermodynamic equilibrium the 
calculations lead to an exact balance between pair production and pair 
annihilation and between Compton and Inverse Compton scatterings. This 
balance can only be achieved If the above stimulation and suppression effects 



23 


are properly taken into account. For systems not In equilibrium, grasar 
action Is possible. We have evaluated. In particular, the absorption 
coefficient for equal photon and particle temperatures and positive particle 
chemical potential (u t >0). For this example of population Inversion, the 
total absorption coefficient can become negative due to the much larger 
probability for stimulated annihilation than for Compton scattering and pair 
production. This type of grasar action can produce a narrow emission line 
peaked at an energy of about 0.43 MeV. This energy Is lower than the peak of 
the spontaneous annihilation emlsslvlty, which occurs at energies greater than 
0.511 MeV. This Is caused by the enhancement of all absorption and stimulated 
emission effects with decreasing energy by the smaller amount of available 
phase space which leads to a larger occupation number for the same photon 
density. In a bath of blackbody photons (u^ ■ 0) and e*-e~ pairs of 
temperature equal to the photon temperature, the threshold for grasar action 
is at w* ■ 1 MeV corresponding to pair densities ~ 10 30 cm“ 3 for T ■ 10®K. A 
temperature of ~ 5x1 O^K is needed to produce this density In equilibrium with 
blackbody photons. 

We have applied our results to gamma ray bursts. In particular to the 
March 5, 1979 transient from which an emission line at ~ 0.43 MeV was observed 
(Mazets et. al. 1979). Similar emission lines have been seen from several 
other bursts as well (Teegarden and Cline 1980, Mazets et.al. 1981). While 
these lines could be produced In a cool skin layer of the source region 
(Ramaty et al. 1980, 1981), grasar action has the advantage of being capable 
of producing a narrow line from a hot and optically thick source and from a 
source region of relatively small emitting area and short duration of line 
formation. But If grasar action Is responsible for the observed 0.4 to 0.46 
MeV emission line seen from gamma ray bursts, then their sources cannot be a 
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neutron star of mass much larger than about 0.6 M 0 . At the surfaces of such 
neutron stars, the gravitational field would shift the line to an energy lower 
than observed. 

There are several difficulties and shortcomings In our treatment. We 
have not shown how the Inversion (it±>0) Is produced. It could In principle 
result from the cooling of the pairs that Is faster than their annihilation, 
or by a rapid external supply of pairs without heating. Cooling by 
synchrotron emission has already been proposed for gamma ray burst sources 
(Ramaty et al. 1980, 1981), but for the high densities that we consider here, 
the required field (8>10^gauss) seems to lead to synchrotron self-absorption 
that could quench the grastr action. We have Ignored other effects of a 
strong magnetic field as v'.ll, by limiting our calculations to Isotropic 
distributions and by using plane wave functions Instead of Landau functions. 
This Isotropic treatment also does not allow the study of beaming effects 
which should be present In a gamma ray maser. Finally, we have not made any 
attempts to study the spatial and temporal development of a system exhibiting 
grasar action. We expect this development to be highly nonlinear. 

We have Indicated In the analytic part of the paper that Compton maser 
action Is possible If the photon temperature Is higher than the electron 
temperature and Indeed In cases where the photons cannot be characterized by a 
single temperature. Such maser action could lead to very Interesting effects 
In the 10 to 100 keV region which we have not yet fully explored. 

As already Indicated, ganma ray burst sources are possible astrophyslcal 
sites where grasar action could occur. The most obvious observational test 
for this would be the observation of a narrow (FWHM « 0.1 MeV) emission line 


at ~ 0.43 MeV. 
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Figure Captions 

1. Density of photons and pairs vs. temperature. The line labelled y 
represents blackbotfy photons. The curves labelled by values of the pair 
chemical potential ^ represent pairs In a Fermi -Dirac distribution. 

2. Emlssivitles and absorption coefficients vs. photon energy In a system In 
thermodynamic equilibrium at 3x1 0®K. The Compton emissivlty (not shown) 
is the difference between the total emissivlty and the annihilation 
emissivlty j^. The Compton absorption coefficient (not shown) is the 
difference between the total absorption coefficient and the pair 
production absorption coefficient K^. The photon and pair densities in 
these conditions are 5.5xl0^^cm"^ and 2.4xl0^cm“3 respectively. 

3. Emissivities and absorption coefficients vs. photon energy in a system at 
3xl(rK having equilibrium between pair annihilation and pair production by 
non-blackbody photons with chemical potential -2.0 MeV. The curves have 
the same significance as in Figure 2. The photon and pair densities in 
these conditions are 2.0xl0^cm”3 and 1 .OxlO^cnT^ respectively. 

4. Emissivities and absorption coefficients vs. photon energy in a system at 
3xl0®K with blackbody photons and an inverted pair population described by 
pair chemical potential 1.1 MeV. The curves have the same significance as 
in Figure 2 except that the Compton absorption coefficient K c is shown 
explicitly. Negative values of total and pair production absorption 
coefficients are represented by dashed curves. The photon and pair 
densities in these conditions are 5.5x10^ cnf^ and 7.3xl0^cm - ^ 
respectively. 

5. Emissivities and absorption coefficients vs. photon energy in a system at 
zero temperature with fully degenerate pairs described by chemical 
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potential 0.85 MeV. The curves have the same significance as In Figure 
4. The photon and pair densities in these conditions are zero and 
1.4xl0 30 cm" 3 respectively. 

6. The development of the intensity of the annihilation line with increasing 
thickness of source. The system is the same as that for Figure 4. The 
labels on successive maxima indicate the thickness involved, and the peak 
energy and FWHM of the line. 
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